Journal of Evolutionary Biology

Eusocial Hymenoptera present the highest known recombination rates among metazoans, which evolved several times independently among bees, ants and wasps. Several hypotheses have been proposed to explain this observation, including stronger selection for recombination caused by coevolving parasites and pathogens, and strong sexual selection among haploid males due to male-biased sex ratios among reproductive individuals. In this article, we explore the effects of haplodiploidy and differential selection between sexes on the evolution of recombination, by analyzing a three-locus model in which selection for recombination stems from negative epistasis between selected loci. Our analytical predictions are compared with the results of individual-based simulations in which deleterious mutations occur along a linear chromosome. Our results show that, at mutation-selection balance for deleterious alleles, increasing the strength of selection against deleterious alleles (due to the effect of male haploidy and/or sexual selection) tends to reduce selection for recombination. However, an increase in the overall magnitude of negative epistasis (which may also be due to male haploidy and/or sexual selection) combined with the fact that recombination only occurs in females may increase selection for recombination substantially. Our model also shows that, in conditions favoring recombination, increasing recombination in meioses leading to parthenogenetic ovules (and male offspring) may yield stronger benefits than in meioses leading to fertilized ovules (and female offspring).

The number of chromosome sets per nucleus is a fundamental trait, but why this number is nearly always two for multicellular eukaryotes remains unclear. Chromosomes are made of nucleic acids, which possess abundant phosphorus (P). Therefore, producing new chromosomes, as well as generating new cells and organismal growth, demands substantial phosphorus. Yet, because P is often limiting in nature, P availability could influence the prevalence of diploidy versus polyploidy. Here, we compare growth rates of diploid and triploid Potamopyrgus antipodarum, a freshwater snail, relative to P availability. Because diploid P. antipodarum are obligately sexual while obligately asexual individuals are polyploid, costs associated with sensitivity to P limitation in polyploids could also help explain the maintenance of sexual P. antipodarum. We raised juvenile diploid and triploid snails on either P-adequate or P-deficient diets and found that independent of P availability, juvenile triploid asexual snails grew faster and harbored higher P content as adults than sexual diploid conspecifics. Together, these results suggest life-history advantages of polyploidy or asexual reproduction that exacerbate rather than ameliorate the cost of sex. These outcomes suggest that P availability is unlikely to be a main driver of ploidy polymorphism or the maintenance of sex in P. antipodarum.

Body size, through its links to various physiological traits, has often been hypothesized to influence evolutionary rates. Negative body size-rate correlations have been reported in the morphological or molecular evolution of several extant vertebrate groups, including mammals, birds, reptiles, and teleost fishes. In this study, we estimated body masses for 89 species of plesiosaurs, a clade of Mesozoic aquatic reptiles, and found that their body size evolution conforms to a three-regime Ornstein-Uhlenbeck process, indicative of constrained evolution. Rates of morphological evolution, inferred using the skyline fossilized birth-death process and the variable-rates model, show minimal support for a correlation with body size in this clade. Our results thus serve as a counterexample, suggesting that the negative body size-rate relationship is not a universal vertebrate pattern, but rather a trend restricted to certain lineages.

Empirical studies have demonstrated that mutation rates may change with individual condition, such as in the case of stress-induced mutagenesis. This has led to the hypothesis that condition-dependent (or "plastic") mutation rates could be selectively favoured, as the increased production of new mutants in times of maladaptation enhances evolvability, the ability to undergo adaptive evolution. However, while empirical evidence for condition-dependent mutation rates is accumulating, theoretical models studying their evolution are lacking. Here, we employ an individual-based simulation approach to examine the evolution of condition-dependent mutation rates in a changing environment. We find that condition-dependent mutation rates consistently evolve when the environment changes at an intermediate pace. Furthermore, populations with condition-dependent mutation rates are substantially better adapted to their (changing) environment. Finally, the evolutionary dynamics of condition-dependent mutation rates are both accelerated and destabilised when the mutation rate is self-referential (i.e., when mutator loci affect their own mutation rate). We conclude that condition-dependent mutation rates (and thus evolvability) can readily evolve in changing environments. Significance statementMutation provides the raw material for evolution. Mutation rates thus tune evolvability, the ability to undergo adaptive evolution: if mutation rates are too low, evolution is impeded; if mutation rates are too high, adaptive traits cannot be maintained. Using a theoretical model, we explore the evolution of plastic mutation rates that systematically depend on the condition of the organism and its environment. An example is stress-induced mutagenesis in bacteria, which is implicated in the evolution of antibiotic resistance. We show that plastic mutation rates readily evolve, providing "well-timed" variation specifically when organisms are poorly adapted. Such plastic mutation rates thus facilitate better adaptation to changing environments, and their evolution provides an example of evolvability itself evolving.

Brain size and structure are hypothesized to be adaptively designed to satisfy the behavioral requirements of securing food and living socially. The importance of these socioecological and sociobiological selective forces in brain evolution is constantly debated. Socioecological divergence is striking in the Neotropical ant genus Neoponera: N. apicalis is a generalist solitary predator forming small colonies of ~100 whereas N. commutata colonies are approximately 10 times as large and workers pheromonally organize cooperatively raids only on Syntermes termite colonies. We interspecifically compared the size and structure of the compound eyes, size and number of antennal glomeruli, mosaic brain scaling and synaptic processing (microglomeruli-MG). Our results indicate that N. apicalis workers have a larger number of ommatidia, antennal lobe glomeruli, and allometrically larger antennal and optic lobes than N. commutata. These sensory traits were associated with differences in higher-order processing architectures in the mushroom body (MB) microglomeruli (MG). N. commutata workers had an allometrically larger MB, perhaps due to their socially complex chemical foraging communication, although MG density in N. apicalis was higher in both the MB lip and collar, regions associated with processing olfactory and visual information, respectively. The increase in MG density in N. apicalis may be associated with higher demands for navigation, learning, and memory, as well as a higher density of antennal lobe glomeruli to support prey odor discrimination. In contrast, N. commutata workers had larger ommatidia and antennal lobe glomeruli. Larger ommatidia correlate with their diurnal/nocturnal habits and a larger MB Our findings indicate that differences in behavioral performance demands associated with socioecological differentiation are reflected in variation in visual and olfactory system structure, brain size, mosaicism, and synaptic organization. Our results support both social and ecological brain hypothesis as drivers of mosaic brain evolution.

Many proposed mechanisms for the evolution of cooperation among unrelated individuals rely on relatively demanding cognitive abilities that are not widespread across taxa. In contrast, individual heterogeneity is a pervasive feature of animal groups, encompassing differences in personality as well as physical and cognitive traits. Such heterogeneity can promote the evolution of cooperation, yet its role has received comparatively little attention, particularly as a source of variation giving rise to social organization such as leadership. A specific form of leadership can emerge under unstable environmental conditions, when some individuals become better suited than others to initiate action and influence the behavior of their peers. Unlike fixed dominance hierarchies, emergent leadership can rapidly adjust to changing environmental conditions, thereby reshaping group organization. Because it does not require the maintenance of stable hierarchies, this form of leadership can arise even in species that do not have the cognitive capabilities to sustain complex social structures. In this work, we investigate the combined effects of individual heterogeneity and emergent leadership on the evolution of cooperation using an evolutionary game-theoretic model in which individuals may assume the roles of leaders or followers according to their strength, representing individual differences in suitability to prevailing environmental conditions. We examine different levels of population heterogeneity together with increasingly complex strategy sets requiring progressively greater informational requirements, allowing individuals to condition cooperation on their own strength, leadership role, or both. Our results show that the interplay between leadership and heterogeneity promotes the evolution of cooperation, particularly when only a small fraction of individuals act as leaders. Under these circumstances, cooperation evolves even when individuals employ the simplest possible strategies. Under harsher ecological conditions, cooperation can be sustained by more sophisticated strategies, specifically by conditional strategies that prescribe cooperation when individuals are strong or leading and defect when acting independently. Author summaryIn this study, we propose that emergent leadership mediated by individual diversity can boost the evolution of cooperation in animal groups. Building on growing evidence on the heterogeneity of animal capabilities and personalities, we focus on the fleeting leadership that emerges in animal groups when facing rapidly changing environmental conditions. We suggest that this type of leadership that emerges from individual differences in strength--a generic quality encompassing those characteristics that make an individual more fit to lead in a given situation--does not require complex cognitive capabilities from the animals and represents a valid alternative to more demanding strategies proposed in the past to explain the evolution of cooperation. Using an evolutionary game theory model, we show that if a population includes a few strong players, these can become influential leaders and guide the actions of their peers to achieve cooperation. Although the naive strategy of always cooperating is sufficient for cooperation to evolve, the introduction of more complex strategies leads players to cooperate only when they are more likely to be recognized as influential leaders. These strategies are more effective in promoting cooperation under unfavorable ecological conditions and are also more robust against exploitation by defectors.

If an individuals niche is determined by its genotype, then competition for limiting resources should be most intense among individuals of the same genotype. Theory predicts this will act to maintain genetic variation, but whether this mechanism operates under natural conditions is unclear. Using long-term observations of a population of free-living Trinidadian guppies, we asked (i) whether competition was strongest between kin, and (ii) whether this process maintained genetic variation. Competition between kin was 1.5-1.8 times stronger than that between non-kin. This contributed to balancing selection: after [~]10 generations, variation was 29% higher than expected under drift. Our results show that relatedness can play a major role in structuring ecological competition, with broader consequences for theories of inclusive fitness. One-sentence summaryHeritable variation is maintained due to resource competition being more intense among kin.

Does the distribution of fitness effects of new mutations vary across the genome? Under the classical Fisher Geometric Model (FGM) we might not expect it to. In FGM, phenotypic traits are envisioned as dimensions of a landscape, with fitness determined by position in the landscape, i.e., the particular combination of traits of an individual. New mutations are represented by vectors that move from an ancestral to a new phenotype. In classical FGM these vectors affect all trait dimensions simultaneously (universal pleiotropy). However, introducing partial and modular pleiotropy into an FGM framework leads to an expectation that parameters of the DFE will vary with mutational pleiotropy-the number of traits affected by individual mutations. Here we address this prediction by investigating whether traits related to mutational pleiotropy, expression level and network connectivity, affect the parameters of the DFE using whole genome data from A. thaliana and C. grandiflora, two closely related Brassica species that vary significantly in their demography and mating system, and therefore, in effective population size and the effects of linked selection. Results were similar across both species. We found that expression level and network connectivity were predictive of the parameters of the deleterious DFE, even once co-correlations among genome biology traits were accounted for. Our results suggest that, across the genome, molecular evolutio(high mutational pleiotropy). nary patterns agree with the predictions of FGM, albeit relaxing the assumption of universal pleiotropy, and that variation in mutational pleiotropy among genes is sufficient to have detectible effects on the DFE. Significance statementHow do the effects of new mutations vary across the genome? If mutations in some genes affect many traits (high mutational pleiotropy), we hypothesise they will be more strongly deleterious, with lower variance in their selective effects. We test this by investigating the distribution of effects of new mutations across genes that vary in features that are related to mutational pleiotropy: expression level, gene network connectivity, and number of associated GO terms. The mean strength and coefficient of variation of selection of new mutations varied across genes with different features in the manner expected by our hypothesis. This demonstrates that important parameters of molecular evolution can vary across the genome with genome architecture.

While asexual species can often outcompete their sexual counterparts over ecological timescales, their long-term evolutionary success is hindered by a diminished ability to purge deleterious mutations and to adapt to changing environments. However, some asexual species persist for millions of years, and a major question in evolutionary biology is how they do so. One solution is to occasionally reproduce sexually, as has been shown in a handful of primarily asexual species. Here, we investigate the possibility of rare sex in the clonal raider ant, Ooceraea biroi. We report the whole-genome sequence of a previously uncharacterized clonal line and, using population genetic and phylogenetic analyses, show that it originated through sexual reproduction between two extensively studied clonal lines. The mitochondrial genome of this clonal line differs from that of the maternal clonal line at only a single nucleotide, suggesting that the sexual reproduction event occurred within the past few hundred years. These results demonstrate that sex occurs sporadically in the clonal raider ant, allowing it to generate new genetic combinations and potentially to overcome some of the costs of asexuality.

High investment in current reproduction can limit future reproductive opportunities, but how selection shapes these hormone-mediated traits remains poorly understood. Androgens can mediate male reproductive investment, and in three-spined stickleback (Gasterosteus aculeatus), exert antagonistic effects on breeding effort versus spermatogenesis. To understand how shifts in reproductive strategy shape this tradeoff, we compared testes transcriptomes and androgen production between two recently diverged stickleback ecotypes that differ in reproductive strategy: the ancestral "common" ecotype, which provides paternal care, and the non-parental "white" ecotype, which has lost paternal care and prioritizes mating effort. During typical breeding, testes gene expression differed little between ecotypes. However, under prolonged summer-like conditions, testes gene expression diverged substantially. Common-biased genes were enriched for meiotic functions and spermatogenic cell type markers, suggesting commons had initiated spermatogenesis while whites had not. Instead, whites expressed higher levels of steroidogenic candidate genes and released significantly more 11-ketotestosterone than commons, indicating sustained investment in current reproduction. F1 hybrids released 11-ketosterone at intermediate rates, suggesting a genetic basis for this divergence. Sustained androgen production in whites may possibly delay the transition into spermatogenesis, limiting investment in future reproduction. These results illustrate how selection on hormonally-integrated traits can drive rapid divergence in life history strategy.

Climate change can threaten population viability by disrupting sex ratios in species whose sex is influenced by temperature. While species with sex chromosomes were historically considered immune, in some species, temperatures can override genetic sex determination via sex reversal, leaving them vulnerable to climate-driven sex ratio shifts. The Tasmanian spotted snow skink (Carinacincus ocellatus), a viviparous reptile with an XX/XY system, provides a compelling case study. While laboratory studies demonstrated that extreme thermal conditions induce female-to-male sex reversal (XX males), its occurrence in the wild remains unexplored, limiting our understanding of actual climate impacts. Integrating 23 years of phenotypic and genetic sexing data across two climatically distinct populations, we provide the first evidence of sex reversal in a wild viviparous reptile. XX reversal occurred in both populations, affecting up to 23.5% of XX births in the warmer population, and was associated with colder minimum daily temperatures. Despite high birth rates in some years, sex-reversed adults were rare. We also identified putative XY females, suggesting bidirectional sex reversal and reinforcing the extreme plasticity of reptilian sex determination. Ultimately, sex reversal could act as an evolutionary trap, potentially compromising population viability as climate instability increases.

Bird migration evolved as an adaptation to seasonally changing habitats. Migratory behaviour can vary within the same species in case of partial migratory behaviour, i.e. one population (or individual) is migratory and another one is resident. Species that exhibit a wide variety of migratory phenotypes provide valuable systems to understand the evolutionary drivers behind different phenotypes and how populations adapt to habitats with distinct seasonality. The European robin (Erithacus rubecula) expresses migratory behaviour in central and northern areas of the species distribution range, whereas populations in the South and on the Macaronesian islands are predominantly resident, providing a suitable system to investigate these questions. We use high coverage whole genome re-sequencing data of 125 European robins to investigate how migration behaviour affects population structure and demography, and how it affects the selection landscape in the genome. Genetic structure in European robins coincides with migratory phenotype and geography and populations are characterised by distinct demographic histories. Our results suggest that both the continental resident population as well as the Macaronesian island populations have derived independently from an ancestral migratory population. Unexpectedly, tests for differential selection revealed extensive positive selection pressure acting across all chromosomes in the resident populations, while selective sweeps are largely absent from migrants. We speculate that this might be an analytical artifact due to mismatching timescales between what population genomics methods can detect and the scale on which migration behaviour likely evolved in the robin. We suggest that future studies on the genomics of migration should more focally account for different time scales on which these processes happen, such as including the wider phylogenomic background of the target species, to capture the full evolutionary history of migratory traits.

In search of recurrent patterns of postzygotic hybrid incompatibility, we investigated the Holarctic 'Old World Swallowtail' (Papilio machaon Linnaeus) butterfly complex, whose members can easily be crossed in the laboratory. How many species this model system comprises remains unclear, as taxa with highly distinctive larvae but uncertain status come into contact with authentic P. machaon subspecies in Southern parts of the Palearctic region. By determining mitochondrial and ITS2 haplotypes within and away from contact zones, we found that these neighboring populations do exchange genes, as expected from F1 hybrids being generally fertile in our laboratory crosses. Nevertheless, recurrent instances of dysregulation of diapause were uncovered in hybrid progenies. In keeping with Haldane's Rule, pupae of the heterogametic (female) sex were either unable to enter diapause or, in reciprocal crosses, unable to resume development after having initiated diapause, whereas F1 males experienced normal, photoperiod-regulated diapause, but passed on abnormal diapause regulation to part of their female offspring when backcrossed. Comparing male and female pupal weights in hybrid progenies provides estimates of parental incompatibility that allow to rank taxa and predict quantitatively the outcome of additional crosses, as expected if the same regulatory system were involved. A survey of the entomological literature confirms that diapausing pupae that cannot resume development ('perpetual nymphs') are a recurrent feature of interspecific lepidopteran crosses. Moreover, of the two parent species of a perpetual nymph, the paternal one generally has fewer broods per year. These observations are discussed in the light of models of the speciation process.

The potential for environmental change to compound the detrimental effects of inbreeding depression in small and isolated populations is a significant concern in conservation biology. Previous evidence syntheses suggested that environmental stress exacerbates inbreeding depression, but were based on limited data. Here, we comprehensively test the relationship between inbreeding depression and environmental stress in natural populations using Bayesian mixed-effects meta-analysis on a large, high-quality data set of 2127 inbreeding depression effect sizes from animals and plants. Our results show that inbreeding depression is significantly higher in benign than in stressful environments. Analyses of both inbreeding depression and stress-induced changes in genetic load supported a unimodal (humped) relationship between the costs of inbreeding and stress intensity, with a peak at intermediate stress. At the highest levels of stress there was, on average, a significantly greater inbreeding load in benign than in stressful environments. We suggest that the lower cost of inbreeding associated with extreme stress results from constraints on the expression of inbreeding depression as fitness and phenotypes decline towards zero. Our findings help to resolve long-standing uncertainty around how inbreeding and environmental change interact, revealing that inbreeding responses vary non-linearly with environmental stress intensity, but showing that stress does not generally amplify inbreeding depression. As such, they will inform both the management of populations of conservation concern and predictions of species responses to global environmental change.

Mutation is the ultimate mechanism that produces genetic novelty, and thus a central ingredient of evolution. Mutation rates are therefore thought to be tuned by natural selection, for example to optimize a delicate balance between the generation of adaptive diversity and the accumulation of deleterious mutations. As this selection occurs over very long time scales, models and simulations have been powerful tools to understand how mutation rate evolves and which factors influence it. Most simulation methods are nevertheless limited by the over-simplicity of the genotype-to-phenotype map they feature, especially regarding the encoding of mutation rate. We modified Aevol, an evolutionary simulator inspired by bacterial genomics with a realistic genome structure and a complex genotype-to-phenotype layer, to allow organisms to evolve genes coding for higher replication fidelity. This setup permits several degrees of realism absent in other models: mutation-rate modifier genes themselves experience a realistic distribution of effects of mutations and diminishing- returns epistasis, similarly to fitness modifiers. Moreover, a lower mutation rate comes with the trade-off of a larger genome to encode the genes improving replication fidelity. We use this setup to test hypotheses regarding the evolution of prokaryotic mutation rate, and its link with genome size and genetic drift. We found that evolution systematically increases replication fidelity, even when this results in lower fitness. We highlight two factors which limit the mutation rate decrease: genetic drift and the supply of gain-of-fidelity mutations.

Hosts defend themselves against parasites through resistance (reducing parasite burden) and tolerance (reducing the fitness cost of infection without affecting parasites). This distinction has important evolutionary implications: resistance is predicted to maintain polymorphism while tolerance tends to fix, and only resistance is expected to provoke parasite counter-adaptation. The reaction-norm framework, which infers tolerance from the slope of host fitness regressed on parasite burden, assumes that a shallow slope reflects parasite-independent host protection. We test this assumption using a within-host model in two variants: microparasites (Model 1, with within-host replication) and macroparasites (Model 2, without). Sublethal immunity impairs the host-exploitation rate of the parasite, reducing both growth and per-parasite virulence without killing them. We show that this generates systematic slope differences among host genotypes that the framework interprets as variation in tolerance. Furthermore, the ranking of slopes across genotypes reverses between linear and sigmoidal damage functions: under linear damage, the strongest immune responder appears most tolerant; under sigmoidal damage, the weakest responder does. Decomposition of the damage reduction shows that virulence reduction accounts for the majority of the effect across both model variants. Thus, the reaction-norm slope cannot determine whether host fitness is maintained by parasite-independent tissue protection or by sublethal impairment of parasites.

Introgression following hybridization is increasingly recognized as a major driver of evolution. However, its importance depends on its frequency in nature, which remains to be quantified. To address this, we provide a snapshot of ongoing introgression in a whole species assemblage (4126 ant colonies). 23% of all 82 local species show signs of introgression, which is more than twice previous estimates. Introgression is typically subtle, yet we find that it contributes measurably to genetic diversity. Species divergence, rather than classical prezygotic reproductive barriers (mating phenology, ecological niche, fine-scale habitat use) constrains introgression, suggesting that the main reproductive barriers are postzygotic at this stage of divergence. Our results indicate that introgression may be a common but often overlooked feature of natural communities.

In a wide variety of animals, developmental crowding results in adults with smaller bodies. The crowding effect on body size in Drosophila melanogaster is canonically attributed to heightened competition for nutrition. However, whether other consequences of crowding also contribute to its effect on size remains an open question. We tested the relative contributions of nutritional competition, oxygen availability, and larval-generated metabolites to the crowding effect on size. We found that while nutrition explains most of the variation in body size due to crowding, oxygen also contributes in a sex- and nutrition-dependent manner. We found no evidence that larval-generated chemicals affect body size. These data confirm a widely suspected but untested role of nutrition in producing the crowding effect on size in D. melanogaster, while revealing an unexpected role of oxygen, and raise the possibility that behavior may be a mediator of density-dependent plasticity. Research HighlightsWe found that both nutrition and oxygen mediate the crowding effect on size in Drosophila melanogaster.

The general random walk model (GRW) of Hunt (2006) is used to infer directional evolution in mean trait values from sparse fossil data by modeling phenotypic change as the accumulated result of small steps with mean step sizes and step variances. Using simulations and real data cases, Ergon (2026) showed that the step variances can be estimated reasonably well only when the mean trait values have small measurement errors, while for fossil data with realistic measurement errors they appear to be extremely difficult to find, and they are often found to be negative. In the simulations Ergon (2026) assumed that the true phenotypic mean values were known. Here, I essentially repeat these simulations under the assumption that only mean trait values with large measurement errors are known, and based on weighted mean squared error (WMSE) comparisons the conclusion is that weighted least squares (WLS) is a better method than GRW. A second conclusion is that WLS is a better method also in the possibly rare cases with large measurement errors where the GRW parameters are estimated well. The GRW method is simply not flexible enough to handle such cases. A third conclusion is that Akaike Information Criterion (AIC) results for GRW models with large measurement errors relative to the step variance may be overly optimistic.

Understanding how social complexity responds to environmental variation remains a longstanding challenge in evolutionary biology. Here, we investigated the drivers of social complexity using intraspecific social variation across seven locations of the obligatory shell-brooding cichlid Neolamprologus meeli in Lake Tanganyika. We quantified the number of subordinate individuals per female territory and examined the effects of predation risk, shell availability, and their interaction. Social complexity increased with shell availability under high predation risk but showed little association under low predation risk. A field manipulative-experiment further demonstrated that increasing shell availability led to higher juvenile retention, indicating a causal effect of territory quality. In addition, removal of subordinates reduced shell availability, suggesting the feedback between group size and territory maintenance. We also assessed genetic population structure based on nuclear SNPs obtained by MIG-seq and found only weak genetic differentiation among localities, suggesting that the observed social variation is unlikely to simply reflect strong genetic subdivision. Together, these results show that predation risk promotes group living, whereas nesting resource availability constrains its extent. Our study highlights that social complexity emerges from the interaction between macro- and micro-ecological factors, providing a mechanistic understanding of the evolution of social complexity and philopatry.